Variety of cosmic acceleration models from massive F{R) bigravity. 



Shin'ichi Nojiri^'^. Sergei D. Odintsov'^''*'^'^, and Norihito Shirai^ 
^ Department of Physics, Nagoya University, Nagoya ^6^-8602, Japan 
^ Kobayashi-Maskawa Institute for the Origin of Particles and the Universe, Nagoya University, Nagoya 464~8602, Japan 

Consejo Superior de Investigaciones Cientificas, 
ICE/CSIC-IEEC, Campus UAB, Facultat de Ciencies, 
Torre C5-Parell-2a pi, E-08193 Bellaterra (Barcelona) Spain 
* Institucio Catalana de Recerca i Estudis Avangats (ICREA), Barcelona, Spain 
* Dept. Gen. & Theor.Phys., Eurasian National Univ., Astana, Kazakhstan 
^ Tomsk State Pedagogical Univ., Tomsk, Russia 

We study the accelerating cosmology in massive F[R) bigravity via the reconstruction scheme. 
The consistent solution of the FRW equations is presented: it includes Big and Little Rip, 
quintessence, de Sitter and decelerating universes described by the physical g metric while the 
^ ' corresponding solution of the universe described by the reference / metric is also found. It is 

\ demonstrated that in general the cosmological singularities of g metric are manifested as cosmologi- 

cal singularities of the reference / metric. However, there are models where cosmological singularity 
does not occur in the universe described by g metric while it occurs in / universe and vice-versa. 
Alternatively, the structure of singularity may change (Big Rip in / universe becomes Little Rip 
in g universe). The consistent solution describing Big and Little Rip, quintessence, de Sitter and 
' decelerating universes for coinciding g and / metrics is established, indicating the connection with 

the convenient single metric background formulation. 

PACS numbers: 95.36.+X, 98.80.Cq 

I. INTRODUCTION 
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Massive gravity history has started from the work [l[ (for recent review, see Q)- It also has been known that 
massive gravity eventually contains the Boulware-Deser ghost Q and vDVZ discontinuity Q in the limit of m — > 0. 
■ The attempts to screen this discontinuity are based on the Vainstein mechanism Q as it was shown, for instance, in 

Ref. on the example of the DGP model 0- 
. It has been realized recently that non-linear massive gravity (with non-dynamical background metric) may 
' be extended to the ghost-free construction with the dynamical metric [lO| . General structure as well as cosmological 
Cn| i evolution of such massive gravity has been studied in Refs. The convenient description of the theory gives 

' bigravity or bimetric gravity which contains two metrics (symmetric tensor fields). Some cosmological solutions of 
, massive bigravity which describe accelerating universe are known . One of two metrics is called physical metric 

fS| ' while second metric is called reference metric. Despite the convenience, the interpretation of such description is not 
easy: indeed, it looks like we have two different universes described by two different metrics. In reality second metric 
just gives the effective description of exotic matter (massive graviton). 

In Ref. we have proposed massive ghost-free F{R) bigravity which is formulated in terms of the auxiliary 
scalars. The cosmological reconstruction scheme of F(R) bigravity was developed in detail with explicit examples of 
accelerating universe. However, there may appear extra ghost (not of the Boulware-Deser type) in such formulation. 
In addition, it is very interesting to understand the connection between background evolutions of two metrics. 

The present work is devoted to the study of accelerating cosmology in massive F(R) bigravity and its relation 
with observable universe. After brief introduction to massive F{R) gravity in next section, we develop the consistent 
reconstruction scheme where extra ghost is not generated. We also study the relation between background cosmologies 
induced by two metrics. We introduce four kinds of metrics, g^jy, gj^^^, /^jy, and f^^. The physical observable metric 
gji^^ is the metric in the Jordan frame. The metric g^i/ corresponds to the metric in the Einstein frame in the standard 
F{R) gravity and therefore the metric g^j^ is not physical metric. In the bigravity theories, we have to introduce 
another reference metrics or symmetric tensor /^j^ and /^^. The metric /^j^ is the metric corresponding to the 
Einstein frame with respect to the curvature given by the metric /^^ . On the other hand, the metric f'^^ is the metric 
corresponding to the Jordan frame. 

Variety of accelerating cosmologies which include Big Rip, de Sitter, quintessence, and Little Rip universes are 
constructed for the space-time described by the metric gj^^ in the situation when the Einstein frame metric is 
fixed. These cosmologies are in good correspondence with observational data which can be shown in the analogy with 
Ref. It is demonstrated that in general, cosmological singularity in physical universe given by the metric is 
manifested in the universe given by the reference metric /^^ or /^^^ and vice-versa (the corresponding observation for 
R massive gravity with matter is presented in Ref. p^). However, we show that there are models where cosmological 
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singularity does not occur in the universe described by the metric 5^^, although it occurs in the universe described 
by the metric /^i, or f^^^. Alternatively, the structure of singularity may qualitatively change, for instance. Little Rip 
in the universe given by the metric manifests as Big Rip in the universe given by the metric f^,y or /^^. It is 
presented the example where the universe given by the metric gjj^^ accelerates while the universe given by the metric 
/aii^ or decelerates. In this sense, the space-time given by the metric /^^ or /^^ (i.e. massive graviton effect) plays 
a role of dark energy in our universe, where the metric is given by gj^^^,. We also propose the qualitative scenario of 
dissolution of the background metric /^^ by the physical metric. The explicit examples of Big Rip, quintessence, de 
Sitter, Little Rip or decelerating universes as solutions of FRW equations in F{R) bigravity when / metric coincides 
with g metric are found by using reconstruction scheme. Finally, general arguments for future singularity occurrence 
in the universe given by the metric f^^ but its avoidance in the physical universe given by 5^^, are presented. 



(1) 



II. F{R) BIGRAVITY 

In this section, we review the construction of ghost-frec F{R) bigravity, following Ref. (for recent review of 
convenient F{R) gravity, see [igI [TtI') . The consistent model of bimetric gravity, which includes two metric tensors 
g^^ and /^j^, was proposed in Ref. [lO[. It contains the massless spin- two field, corresponding to graviton, and massive 
spin-two field. The gravity model which only contains the massive spin-two field is called massive gravity but we 
consider the model including both of massless and massive spin two field, which is called bigravity. It has been shown 
that the Boulware-Deser ghost Q does not appear in such a theory. 

The starting action is given by 

S-bi =Ml j d-^xy^-dctgR^^^ +M] J d^x^^- dct f R^^^ 
+ 2m^M^g / d^x^-detgJ2Pn e„ 

n=0 

Here i?*^^^ is the scalar curvature for 17^1, and R^^'' is the scalar curvature for /^j^. Afoff is defined by 

Furthermore, tensor \/ g^^ j is defined by the square root of g'^'^fpi^, that is, (^\/g^^f^ (^\/9^^f^ ~ 9^'' fpv- For 
general tensor A^^,, e„(A)'s are defined by 

eo(A) = l, ei(A) = [A], e2(A) = i([A]2 - [A^]) , 
e3(A) = i([A]3-3[A][A2]+2[A3]), 

e,{X) = ^([A]4 - 6[A]2[a2] + 3[AT + 8[A][A3] _ Q[X^]) , 

ek{X) = for fc> 4. (3) 

Here [A] expresses the trace of arbitrary tensor X^^^: [A] = A'^. 

In order to construct consistent F{R) bigravity, we add the following terms to the action ([T|): 

S^^- A/2 I d^x^/^d^ I ^g'^-'d^^d,^ + V{ip) I + y" d^xLn.atter ie'^9^., '^^) , (4) 

Se = - M] j d^x^/^d^f ^rd^id,^ + U{0 } • (5) 

By the conformal transformations g^^, — > e^^^g'J^^, and /^^ c^^/^^^, the total action Sp = Shi + S^ + S^ is transformed 
as 



Sp =M] I d^x^^dctp {e--'/?J(-^) +e-2«C/(0} 



. 4 

2m2A/2ff / d'^xyj- dot /3„c(*-2)'^-7«f 
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The kinetic terms for ip and ^ vanish. By the variations with respect to ip and ^ as in the case of convenient F{R) 
gravity [ill, '^^ obtain 

=2m'M!f, E (f - 2) e(*-^)'^-^«e„ Ug-^-'fA + {-c-^i?^^^) - 2e-^^V{^) + o-^^V'{^)] , (7) 



The Eqs. ([7]) and ([8|) can be solved algebraicahy with respect to and as (p = f i i?'"^^-', i?'"--^-', e„ ( y g'' ^.f-' 



and £, = ^ \^R'^^^\R^^^\en 9"^ ) )' ^^^^^^^uting above ip and ^ into dH), one gets bigravity: 



(9) 



(10) 

Note that it is difficuh to solve Eqs. ([7|) and ([8|) with respect to ip and ^ explicitly. Therefore, it might be easier to 
define the model in terms of the auxiliary scalars tp and ^ as in 

III. COSMOLOGICAL RECONSTRUCTION 

Let us consider the cosmological reconstruction program following Ref. [l^ | but in slightly extended form. 
For simplicity, wc start from the minimal case 

^bi ^^^g J d'^x^J-detgR^s^ + M] j d*xy/~ dot f R'-^^ 

+ 2m'^M^ff I d'^xy/- det g (s - tr y/g-^f + det y/g-\f^ . (11) 



In order to evaluate Sy^g ^f, two matrices M and N, which satisfy the relation — N arc taken. Since 6 MM + 
MSM SN, one finds 

ti-SM = itr (M-^SN) . (12) 
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For a while, we consider the Einstein frame action (|TT|) with (|4]) and ^ but matter contribution is neglected. Then 
by the variation over g^^, we obtain 



+ Af 



#2 



1/3 \ 3 ' 



(13) 



On the other hand, by the variation over /^^, we get 



■Mi 



(14) 



We now assume the FRW universes for the metrics g^jy and /^i/ and use the conformal time t for the universe with 
metric 5^^. ^ : 

3 / 3 \ 3 3 

dsl = gf^^dx^dx" = a(tf -df + ^ {dx'f , dsj = ^ ff.^dx^dx'' = -c{tfdf + h{tf ^ {dx'f . (15) 

Then (t, i) component of gives 
and (iii) components give 



(16) 



= Afg2 (^2ij + + TO2M<32ff (Sa^ - 2ah - ac) 



Here iJ = d/a. On the other hand, {t^t) component of p4|) gives 



= -3M'jK^ + m^M^ffC^ [ -3 + — + — 

a a 



and components give 



7b 



= Af? 2X + 3/^2 - 2Lif + mHd^ffC^ 3 + --^^ + -U{Oc{tf M 



a a 



3-9 1, 



(17) 



(18) 



(19) 



Here K = b/b and L = c/c. We now redefine scalars as (p = f{rj) and ^ = ^(C) and identify 77 and C with the conformal 
time t, rj ^ ( ^ t. Hence, one gets 



Here 



w(i)A/| =4Af2 (^i/ - H^j + m^M^ff (-ac + ab) , 
F(t)a(t)2A/2 = - (2H + AH^) - m^M^ff (6a^ - 5ab ~ ac) , 



a{t)M'j =AMj (k - LK) + 2m^MlsC^ 



c b 
a a 



M] (2k + 6K^ - 2LK^ - m^M^ffC^ ( 6 



3c 7b 
a a 



u;{v)^3cp'if^)\ f(r;) = F(^(r?)) , a(C) = 3C'(C)^ C/(C) = C/ (C (C)) 



(20) 
(21) 

(22) 
(23) 

(24) 



Therefore for arbitrary a{t), b{t), and c{t) if we choose uj{t), V{t), cr{t), and [7(i) to satisfy Eqs. (|20ll23p . the cosmo- 
logical model with given a(i), b{t) and c{t) evolution can be reconstructed. 



^ In Ref. O . we have used the cosmological time instead of the conformal time. The use of the conformal time simplifies the formulation. 
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IV. ACCELERATING COSMOLOGICAL MODELS 

Let us construct some examples of cosmological models which describe Big R ip universe [l^, [2^, quintessence 



universe, de Sitter universe, decelerating universe, and Little Rip universe |21|-|23|. The physical metric, where the 
scalar does not directly couple with matter, is given by multiplying the scalar field to the metric in the Einstein frame 
in gl) or dni): 

9i. = e'^Sp. • (25) 

In the bigravity model, there appears another (unphysical) metric tensor f^^, besides g^^. In our model, since the 
matter only couples with g^j^^, the physical metric could be given by ^^j^ in ((25|) . In principle, however, there could 
be a matter coupled with /^i,. Then one may consider two space-times, one is described by 5^1/ (or p^^,) and another 
by //ii/- Although the matter in the space-time described by /^^ could not directly couple with the matter in the 
space-time described by f^^ji^, the matters can interact with each other via the propagation of massless and massive 
gravitons. Furthermore, as clear from Eqs. and the metric /^^ affects the geometry of the space-time 

described by g^^, and vice- versa. 

In this section, we show that in general, there occurs future singularity in the space-time described by g^i, when 
the singularity occurs in the space-time described by f^^. We also find some examples where the future singularity 
does not appear in universe described by 17^1/ (more exactly by gj^^^) even if future singularity appears in universe 



described by This might be a counterexample for the observation of Ref. |24| where R bigravity with matter 
was considered. Although the physical metric gj^^^ in g universe is given by (j25|) . it is not so clear what could be a 
physical metric in / universe since it depends on the coupling with matter. Anyway we may consider both of /^^ in 
the Einstein frame and the metric 

/m. = , (26) 

in the Jordan frame. 



A. Conformal description of accelerating universe 

In our formulation, it is convenient to use the conformal time description. Hence, let us describe how the known 
cosmologies can be expressed by using the conformal time. Especially we present the explicit expressions of the de 
Sitter, phantom, quintessence, decelerating, and also Little Rip universes. 

The conformally flat FRW universe metric is given by 

ds^ = a{tf (-df- + {dx'Y] ■ (27) 



Eq. (Uni) with (HIl) shows 



that is. 



Using ([24]), we flnd 



Here H = i§. 

a at 



e'^(*)a(t)2 = d{tf , (28) 



ip ^ -2 In a{t) +2n In d(t) . (29) 



uj{t) ^U(H - hV . (30) 



In Eq. (P7|. when a(t)^ = j2i the metric (P7)) corresponds to the de Sitter universe which may describe inflation or 

dark energy in the model under consideration. On the other hand if a(i)^ = ^ with n 7^ 1, by redefining the time 
coordinate as 

di=±—dt, (31) 
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that is, 



the metric (l27l) can be rewritten as 



m 

i=± -t^-\ (32) 

n — 1 



ds^ = ~dP + (^±{n ' " ^ {dx'f . (33) 

Eq. p3p shows that if < ri < 1, the metric corresponds to the phantom universe, if n > 1 to the quintessence 
universe, and if n < to decelerating universe. In case of the phantom universe {0 < n < I), one should choose + 
sign in ± of (|3T|) or (|32|) and shift i in ([33)) as f — > f — to- Then i = to corresponds to the Big Rip and the present 
time is t < to and t — > oo corresponds to the infinite past (t — > — oo). In case of the quintessence universe (n > 1), we 
may again choose + sign in ± of ([5T|) or ([5^ . Then t corresponds to t — ^ +oo and t — )■ +oo to t ^ 0, which may 
correspond to the Big Bang. In case of the decelerating universe (n < 0), we may choose — sign in ± of ([51]) or (|32p . 
Then t ^ corresponds to t ^ +oo and t +oo to t ^ 0, which may correspond to the Big Bang, again. We should 
also note that in case of the de Sitter universe (n = 1), t — >■ corresponds to f — > +oo and t — ±oo to f — — oo. 

One may also consider the Little Rip universe, where there is no future singularity but the Hubble rate in terms 
of the cosmological time t becomes infinite when t goes to infinity. As the universe expands, the relative acceleration 
between two points separated by a comoving distance I is given by ^ (1/a) [d^d/dP), where a is the scale factor. An 
observer at comoving distance I away from a mass m will measure an inertial force on the mass of 

a dt'^ 

Let us assume the two particles are bound by a constant force Fo- If -Finer is positive and greater than Fo^ the 
two particles become unbound. This is the "rip" produced by the accelerating expansion. It leads to finite-time 
disintegration of bound objects much before the singularity. 
An example is given by 

i§=goAe^*-A. (35) 
a at 

The last term —A is added for the convenience in the explicit calculation but this term can be neglected for large t 
compared with the first term. By the results of the Supernova Cosmology Project [2^, the parameter A in ([55]) is 
bounded as 

2.37 X IQ-^ Gyr^i < A < 8.37 x 10"^ Gyr^^ . (36) 

Eq. dSS]) gives 

'=-1—^ 5 = e^°°"-^* = - (37) 
Ago Xt\n(Xgot) 

By using p4p . we find the inertial force is given by, 

Finer = ml jsoA'c^* + (goAe^* ~ ^)'} ' 



(38) 



1 da 1 , ^ , , , 

^ = ^-37--7 1 + ■ (39) 



which goes to infinity at infinite t. 
Combining ([55)) and ([57|. one finds 

a dt t \ In (A^o^) 

Thus we get an explicit example of d and H corresponding to the specific Little Rip universe. 

If the space-time described by the metric g'j_^^ describes the universe where we live, the functions c{t) and h{t) are 
not directly related with the expansion of our universe since the functions c{t) and h{t) correspond to the degrees of 
freedom in the Einstein frame metric /^y. 

Therefore one may choose c(t) and b{t) in the consistent way convenient for the calculation. This does not mean 
c{t) and h{t) are not relevant for the physics besides the expansion of our universe since the matter in the universe 
given by the metric /^j^ or J^^,, if any, weakly interacts with the matter in our universe via the massless and massive 
gravitons. In the following, we consider two choices of c{t) and 6(t), that is, a{t) = c{t) = 1 case and a{t) — c{t) = h{t) 
case. Note that the last case probably helps to simplify the formal description of the theory, indeed reference metric 
seems to be dissolved by physical one. 
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B. Dark energy universe with a{i) = c{t) — 1 

In this section, making the choice a(t) = c{t) = 1, we exphcitly construct Big Rip (phantom), quintessence, de 
Sitter, decelerating or Little Rip universes. 

When a{t) = c{t) = 1. the Einstein frame metric g^,^ expresses the flat Minkowski space although the metric we 
observe is given by g'l^. Eqs. (HO]), ([21]), (EH), and ([231) with ^ are simplified as follows, 

UM^H^ = - m^M^ff {1-b) , (40) 

V{t)M^ =bm^M^g (!-&) = eOM^H"^ , (41) 

a{t)Mj =4:Mfk + 2m^M^ff (1 - ^) = iM^K - 24:M^H^ , (42) 

U{t)Mf = - A/^ (2k + 6K^^ - m^Mls (3 - 76) = -Mj (2k + GA'^^ + 'im'^Mlf^ + SAM^H^ . (43) 

Since Eq. (|i(71) can be solved with respect to b as 

2 



wc find 



off 



Then we obtain 



i 2 



7J2 



(45) 



ait)Mf = ^ ^ ^ - 2AM^H^ , (46) 

(1 + ;^^^) 



Umkf = ^ ^^^^^ ^ ^ + 4m^A4^ff + 84M,^i/^ . (47) 



Note that a{t) is not always positive, which conflicts with Eq. (P^ . Therefore the model cannot be identified with 
the analogue of the F{R) gravity if a{t) is not positive. One should also note that negative t7(i) also means the scalar 
field ^ or C is a ghost although there does not appear the Boulware-Deser ghost 



1. Construction of the models describing Big Rip, quintessence, de Sitter and decelerating universes 

As shown in Subsection IIVAI Big Rip, quintessence, de Sitter and decelerating universes are described by the scale 
factor a{t)^ = in terms of the conformal time t. Let us construct the models with the scale factor a(t)2 = j^. 
Studying the properties of such models, we show that there exist models where future singularity does not appear in 
the space-time described by gj^^^ although there appears future singularity in the space-time described by /^^ or f^^. 

Starting from (pS)) . one finds 



2AM'in^ 



a{t)M] = 2« V ,SiL^ ^ 48 

^ ' i ' 12A/>2 s2 ^2 ' ^ I 



1 



liiMfM^n'' /, 28M„n^ 



A4fff V 



2^2 



U{t)Mf. = ^ ' + 4m' M^^ 1 + 21 . . (49) 
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When t is small, a{t) behaves as 

(50) 




In order to avoid the ghost and for the consisteney, we obtain a constraint as follows: 



^>3n^ (51) 
g 



On the other hand, when t is large, the second term dominates in Eq. 

Therefore, a{t) becomes negative and conflicts with ([M|) and the large t region cannot be described by F{R) bigravity 
with the choice a{t) = c{t) = 1. Anyway the model with pS)) is consistent for small t as long as the condition ([?T|) is 
satisfied. One may consider the model (|48p as an example where there occurs singularity in the universe described by 
ffiu or /^j^ but there does not occur any singularity in g universe described by gj^^. Note that the universe described 
by 5^1/ is the fiat Minkowski space-time. For the simplicity, we now only consider n = 1 case, which corresponds to 
the de Sitter space-time in ([?7| for the metric gj^^. The de Sitter space-time describes the acceleratingly expanding 
universe and does not show any singularity. One considers the case t — > 0, which corresponds to infinite future, where 
ghost does not appear and where F{R) bigravity description is valid. Since H = d/a = 1/t, Eq. (|44| shows 



12M2 12M2 



err cr± 



(53) 



Using (|15p . since a{t) = c{t) = 1, we find 



dsj = f^.dx^dx" = -c{ifde + b{tf {dx'f ~ -dt" + / 1 {dx'Y . (54) 

ti.i^=0 i=l \ off 



Furthermore, Eq. (f50| and a{t) in ([24|) indicate 



t 2_ 1 / 8Af2 



e~±2aln-, ^^-12 . (55) 

to 12 \ m-'M^ff I 



Here to is a constant of the integration. Then 'physical' metric f^^ in ([^5)) is given by 

2 1 




{d'^if = E fUd^'dx^ = ^'dsj cx I -dt^ + ( ) (dx^^ | . (56) 



Note that there could be two models corresponding to the signature ±, which give the identical cosmological expansion 
in the space-time described by gj^^. Furthermore we define the cosmological time tj in the space-time described by 

fU by 



l±a 



Then the metric has the following form: 



{dsjy cx -dtj + ' {dxf . (58) 

TO^M^f 1 1 ± a| tj^ 

Hence, if > 0, tj ~ corresponds to Big Rip singularity and if < corresponds to Big Crunch or Big 

Bang singularity. In summary, we presented the model where there occurs cosmological singularity in the universe 
described by f'^^ but there does not occur any singularity in the universe described by gj^^. 
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2. Little Rip universe 



Let us discuss Little Rip universe. In this section we show that the Little Rip cosmology in the space-time described 
by the metri c ff^ ^^ corresponds to the Big Rip universe in the space-time described by the metric /^j,. When t — > 0, 
substituting (|^ into (|3t]), we find 

b{t) ^ . 59 

Then 

K — ^ , k ^ . (60) 

tln{Xgot) t2(in(AgoO) 

Using third equation in and (jH]), one obtains 

3C'(t)2 ^ik 2 ' (61) 

which can be integrated to give 



Hence, the metric /^^ is given by 



C=^/-ln|ln(A5ot)| . (62) 



d^f = E fUd^'dx-^ = ^^ds) ^ |ln {\g,t)\^^ [ -dt' + ( —^^^—^ ^ {dxf ] . (63) 

\P {hi {Xgot}) I .^^ 



fi.v—O 

The cosmological time t/ is defined by 



/■* 

tf= / df|ln(A5o0r^ ■ (64) 



The integration in the r.h.s. is given by the incomplete F function and is finite when t — >■ 0. Therefore t = 
corresponds to the finite tf. When t = the scale factor diverges 

ait)f ^^^1-^. (65) 

i2|ln(Affot|)'"^^ 

Therefore, finite-time future singularity occurs at t = 0. Hence for the case a{t) — c{t) — 1, the Little Rip in the 
universe described by corresponds to the Big Rip in the universe described by f^^^. 



C. Dark energy universe with a{t) = b{t) = c{t) 

As we observed above, general bigravity formally describes two gravities which are related via some coupling term. 
Generally speaking, this is rather formal presentation as second universe described by reference metric is kind of 
effective description of exotic matter (massive graviton). Nevertheless, it may be useful to clarify the reference metric 
in better way. To do this we propose that in the course of the evolution the second universe metric may become equal 
to the physical universe metric (of course, perturbation theories are different). In other words, / metric is dissolved 
by g metric. After that the future background evolution is conveniently described by the single metric object. 

Let us choose a{t) = c{t) = b{t) and therefore H = K = L. From ([20]) and ((22|) . we find ijj{t) = a{t) and therefore 
ip{t) = ^{t) although V{t) ^ U{t) from ^ and which tells, not only g^^ = f^^, but gf^^ = f^^ from ^ and 
(|26p. Hence if there is any singularity in the space-time described by f^^y or fj^^, there appears an identical singularity 
in the universe described by (7^^ or gjj^. For the choice a{t) = c{t) = 1 as in Subsection IIVBI there may appear 
ghost in general and the model cannot be interpreted as a model of the F{R) bigravity. However, for the choice 
a{t) = c{t) = b{t) in this Subsection, there does never appear the ghost as it will be shown below. 
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By choosing a{t) = c{t) = b{t), Eqs. ([201), ([21]), ([221), and ([231) are simplified as 

i {h - Hj' =H - H"^ , (66) 

V{t)a{tfMl = - Ml (2H + AH^'j , (67) 

a{t)M] =4M| (^H - H^^ , (68) 

U(t)a{tfM'j = - Mj (2H + AH^'j + Am^M^gO^ . (69) 
Here, pop is used. Comparing with (|68p. we find 

a{t) = uj{t) ^ (^H - hY . (70) 

Since ct(<) is positive, there does not appear ghost, which is different from the case of Ref. [ij] as is explained in 
Subsection IIV CI In order to determine a, however, one should solve the differential equation (|66p . In case that it has 
no solution, we cannot reconstruct such a model. 



1. Construction of Big Rip, quintessence, de Sitter and decelerating universes 

Let us consider the construction of the models, which describe Big Rip, quintessence, de Sitter and decelerating 
universes, where the scale factor is given by a{tY = 
In case a{t)'^ = jtk, the solution of (|66p is given by 



„ ho , -6n-l±V-12n2 + 12n+l ^^^^ 
^ = y , ho = . (71) 

In order for ho to be real, the parameter n should satisfy the following condition: 

3-2V3 3 + 2V3 

^—<n<^—. (72) 

Then when > n > = -0.07735 • • • , the universe is decelerating, when < n < 1, the universe is in phantom 

phase, and 1 < n < ^+^'^ 1.07735 • • • , the universe in quintessence phase. From Eq. ((7T|) it follows 

a = aot''" . (73) 
Here oq is a constant. Then from Eqs. (|66| . ([Ff]). (|68)) . and (|69|. we find 



.,,,.1^, ,«,.1H^1|±^, n„-^. .,o^*^.^^^ m 

Thus, we have constructed the models, which describe Big Rip, quintessence, de Sitter and decelerating universes. 
Moreover, / metric seems to be dissolved by g metric for such the background evolution. Note that we can present 
the dynamical solution with similar properties. Starting from / universe where b{t) = c{t) one can reconstruct the 
evolution where these parameters just approach to a{t) so that they happen to coincide from some specific moment 
(for instance, pre- inflationary stage). Such a moment will correspond to dynamical background dissolution of / metric 
(the perturbation theories are still different). We will not present here the corresponding results because they are 
quite complicated. 



2. Little Rip universe 

The next example is Little Rip universe unifying again the description of both of cosmologies given by g"^^ and f^^. 
When t ~ 0, Eq. (p9| indicates H ^ ^ tin{\got) ■ Then the solution of Eq. ([66|) is given by 

H = 2+o( 3). (75) 

t(ln(A5oO) \t{\ni\got)f 
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Then since H <^ H when t ^ 0, one gets 

u;(t) = ait) ^ H"^ ^- 7 , V ^il ^ -2H - . (76) 

^' ^' t^hiiXgotf i^(ln(A,9oO) 

Thus, the background cvohition of both of the universes described by and fj^^, happens to coincide as Little Rip 
cosmology. 

D. General arguments for future singularity occurence 

We may consider more general case that there occurs singularity in the universe described by /^^ but there does 
not occur any cosmological singularity in the universe described by gj^^. 

• Even if b{t) and c(t) in the metric (|15p are not singular, if ^ is singular, the metric /^j^ ([26)) becomes singular. 
In this case, pO|. ([2T|) . and p4)) show that a{t) and (/)(i) and therefore p^i/ and gj^^, can be non-singular. 

• Even if h{t) and c{t) have a singularity, if bit) — c{t) is non-singular, (pil)) and (|24p show that a(<) and </j(t) and 
therefore g^v and can be non-singular although (f2T|) indicates that l/(t) is singular. 

Let assume h{t) — c{t) has a singularity at i = 0. Then (|20p and ([^^ demonstrate that a(t) and/or (y5(i) are 
singular at t = 0. If a{t) ([25)1 behaves as 

a(t) ~ ^ , (77) 

the metric gj^^, describes (asymptotic) de Sitter space and therefore g universe is not singular. Using (|77p . 
eq. (PO)) can be written as 

12 (1 + ^ + i7(f)2) =a(h-H^)+ (6(0 - c(i)) a{t) . (78) 



Then one can consider the following cases: 

— When a{t) and therefore Hit) are regular at i = 0, if 



12Af2 / 1 2i7f0) 1 



the universe described by gjj^^ is not singular even if the universe described by /^^^could be singular. 

— Let assume that when t ^ 0, a{t) behaves as a{t) ^ aot'^° with constant oq and Hq. Then Eq. ((78)) 
demonstrates if b{t) — c{t) behaves as 

bit) - c{t) ~ ^ ° ' , (80) 

the space-time described by g^j^ is not singular. 

It looks that all possible cases where the space-time described by f^i, or fj^^ has a singularity but there does not occur 
any cosmological singularity in g universe are presented in this subsection. 

V. DISCUSSION 

In summary, we studied massive F{R) bigravity in the conventional description with two metrics. The variety of 
cosmic acceleration cosmologies is found as explicit solution of FRW equations. In particular. Big and Little Rip, 
de Sitter, quintessence and decelerating universes are constructed for Jordan frame physical metric g when Einstein 
frame metric g is fixed and corresponding solution for reference metric / is also presented. The relation between 
properties of g and / cosmologies is investigated in detail. For instance, it is demonstrated that, in general, the 
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physical g cosmological singularity is manifested as metric / cosmological singularity. However, there are examples 
where cosmological singularity of physical g universe does not occur in the universe described by reference metric / 
and vice-versa. Furthermore, the structure of singularity may qualitatively change: what looks like Big Rip in one 
space-time manifests as Little Rip in other universe. The solution of FRW equations where two metrics just happen 
to coincide is presented for Big and Little Rip, quintessence, de Sitter and decelerating universes. In this case, the 
background evolution is described via single metric which looks quite convenient even keeping in mind that second 
metric is just the effective description of exotic matter (the perturbation theory for two metrics is also different). 

It could be interesting to consider how we can observe the manifestation of the space-time described by the metric 
Itiv or /^j^. As we mentioned, we can observe the matter via massless and massive gravitons. Therefore not only dark 
energy but also dark matter might be a matter in the space-time given by the metric or /^j^. 
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